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Abstract 

This is a study of the path probabiHty by numerical simulation of 
stochastic motion of non dissipative or quasi-Hamiltonian systems. 
This ideal dynamical model implies that, apart from the random 
forces, the system is only subject to conservative forces or that, in 
the presence of friction force, the dissipated energy is negligible with 
respect to the variation (work) of the conservative force. In the nu- 
merical experiment, we used small particles subject to conservative 
forces and to a Gaussian noise (random displacements). The path 
probability was determined by observing a large number of particles 
all moving along many sampled paths between two fixed points in 
a fixed time period. It is found that the path probability decreases 
exponentially with increasing action of the paths. 
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1 Introduction 



The path (trajectory) of stochastic dynamics in mechanics has much richer 
physics content than that of the regular or deterministic motion. A path 
of regular motion always has probability one once it is determined by the 
equation of motion and the boundary condition, while a random motion may 
have many possible paths under the same condition, as can be easily verified 
with any stochastic process including Brownian motion[l, 2, 3, 4]. For a 
given process between two given states (or configuration points with given 
durations), each of those potential paths has some chance (probability) to 
be taken by the motion. The path probability is a very important quantity 
for the understanding and the characterization of random dynamics because 
it contains all the information about the physics: the characteristics of the 
stochasticity, the degree of randomness, the dynamical uncertainty, the equa- 
tions of motion and so forth. Up to now, only a few theoretical works have 
been devoted to the study of this probability. Examples include the appli- 
cation of the Feynman path integral technique to Brownian motion[l] and 
the large deviation theory['>, G] which introduced an pseudo-action or rate 
function to characterize the path probability. To our knowledge, no exper- 
imental work or numerical experiment has been made to measure the path 
probability. This is certainly related to the difficulty of experimental observa- 
tion of a large number of stochastic motions which is necessary to determine 
correctly the path probability. The purpose of the present work is to make 
numerical experiments to study path probability and its relationship with 
the conventional mechanics quantities such as position, velocity, energy and 
action. 

The present paper describes the first results of this work. At this stage 
we focus on the study of the stochastic motion of non dissipative or quasi- 
Hamiltonian system. This means that, apart from the random forces and 
energy fiuctuation, the systems contain only conservative forces without en- 
ergy dissipation. The average energy of the system can or can not change 
during the entire period of the motion. This seems an ideal motion. But real 
systems of this kind do not lack in nature. Any fiuctuating nonequilibrium 
process in an isolated system subject or not to conservative forces (expan- 
sion of isolated gas, heat fiow in isolated solids or liquids and so on) can be 
included. Weakly damped motion can also be included if the energy dissi- 
pated is negligible compared to the variation of potential energy during the 
motion governed by the conservative force and fiuctuation. It is not difficult 
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to create these approximate ideal motions in laboratory. We would like to 
mention that this work has been motivated by a recent theoretical extension 
of Hamiltonian and Lagrangian mechanics to a stochastic formalism[7, 8]. 

It should be stressed that the stochastic motion of those systems which 
are statistically in equilibrium with environment, or whose motion is strongly 
damped such as the usual Brownian particles after long moving time, is not 
considered at this stage. The membership of these systems in the family of 
the Hamiltonian/Lagrangian mechanics, even in the limit case of vanishing 
fluctuation, is still an open question which still nourishes a constant effort [9, 
10, 11]. In order to tackle strongly damped stochastic motion, the model we 
use in this work is to be complemented by the consideration of dissipative 
force and energy. 

The model of our ideal system is a small silica particle subject to conser- 
vative forces and a random force yielded by white Gaussian noise imitating 
the random forces in air or water at room temperature. The Gaussian noise is 
chosen at this stage since the time scale of each step of the simulated motion 
is larger than 10~^ s. The magnitude of the noise can be controlled to ensure 
that the motion is sufficiently smooth and that the velocity defined by the 
displacement at each step is as good as the instantaneous one of the motion. 
In the real random dynamics, the instantaneous velocity has been experimen- 
tally measured with sufficiently small measuring time scale[12, 13, 14]. The 
measured result should be more and more precise with smaller and smaller 
scale. This is certainly an experimental argument for the use of velocity in 
the Langevin equation and the Ornstein-Uhlenbeck model[15]. In very small 
time scales < 10~^ s, the noise is possibly non Gaussian. The work with non 
Gaussian noises will be described in a later report. 

2 Technical details of numerical computation 

We consider a large number (~ 10^) of silica (Si02) particles of mass m = 
1.39 X 10~^^ kg undertaking one-dimensional stochastic motion in conserva- 
tive force field. The spherical particles with l-/im-diameter, moving in air 
or liquid water under Gaussian noise, move from the initial point Xa to the 
final point Xb, over a given n steps, through different paths, i.e., different 
sequences of random positions {xa,Xi,X2 ■ ■ where Xj is the posi- 

tion at time ti with Xa=xo and Xb=Xn- We chose ?7,=10 (due to the limited 
computation time) with equal time increments dt = U — ti-i. Different values 
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of time scale will be used from dt = 10~^ s to dt = 10"^ s. 

The noise or random displacement at tj is given by the following distri- 
bution ^ 

P{xi - Xi_i,U - ti^i) oc e"''^'(*»-'-i> , (1) 

where D is the diffusion constant. If the particles are free, their paths 
{xa, Xi, X2 - ■ ■ Xn-i, x^} Can be generated by a Gaussian noise in the following 
way: 

Xi = Xi_i + axi, (2) 

where a = y2Ddt is a parameter adjusting the magnitude of the random 
displacement Xi — and Xi is a Gaussian-distributed noise generated by 
the computer at each step i. 

For particles in conservative force fields, the model consists in separating 
the motion into two parts: a random part given by Eq. (2) and a part 
described by the solution Ui = f{ti) — of the Newtonian equation of 

motion under a conservative force. The total displacement of each step is 
then given by 

Xi= Xi-i + axi + f{U) - f{U-i)- (3) 

It is obvious that, in the case of vanishing noise, the motion of Hamilto- 
nian/Lagrangian mechanics is recovered. 

The left panel of Figures 1-5 illustrates some sample paths generated 
by Eq. (3) for 5 kinds of conservative forces. The samples are around the 
least action path with sufficiently different actions from the least one. In the 
simulation, each sample path is in fact a bundle or a tube of a small thickness 
5. The larger 5 is, the more particles will go through each path from a to h. 

The numerical experiment consists in observing the total number of par- 
ticles N moving from a to 6 through whichever path and the number of 
particles Nk moving along a given sample path k from a to h. The prob- 
ability that the path k is taken is determined by Pk = N^/N (with large 
N). 

For each sample, we calculate the action (called Lagrangian action 
from now on) and the time integral of Hamiltonian (called Hamiltonian action 
Ah) for the sake of comparison. The velocity at time step i is calculated by 
"^i — xhe kinetic energy is Ki = \mvf, and the two actions by 

Al = Etlllrnvf - V{x,)] ■ dt and Ah = EZA^mvf + V{xi)] ■ dt. The 
magnitudes of the displacements and forces are chosen such that the kinetic 
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and potential energy are of the same order of magnitude. This allows to 
clearly distinguish the two actions along a same path. 

Simulations were performed with 5 potential energies: free particles with 
V{x) = 0, constant force with V{x) = mgx, harmonic force with V{x) = 
i/cx^ and two other higher order potentials V{x) = ^Cx^ and V{x) = \Cx^ 
{C > 0) to check the generality of the results. The result of the simulations 
are presented below. 

3 Determination of path probability distri- 
butions 

Once the path probability is determined as mentioned above, one can search 
for its correlation with the actions (Lagrangian or Hamiltonian one) by draw- 
ing the values of the probability of the sampled paths against the two actions. 
The results are given below for each potential. 

3.1 Free particles 

Free particles has zero potential energy. So these is no difference between 
the Lagrangian and Hamiltonian actions. From the right panel of Figure 1, 
it is obvious that the probability for a path k to occur is given by 

Pk{A) = |e-^^^ (4) 

where Aj. is either the Lagrangian or Hamiltonian action of a path k between 
the a and b. The slope gives 7 ~ 6.7 x 10^^ J~^s~^. The normalization 
function Z can be analytically determined by normalization with the path 
integral technique [16] 

1^ (h^ 1^ dx2 1^ dx^ ^ 1^ cix„_i ^^^^^^^ ^ ^ 
J~oo 6 J~oo 6 J-00 6 J~oo 6 

with fixed Xa and Xb, or numerically by the value of lnp{A = 0) which can 
be found with the distribution curves in the figures. 

3.2 Particles under constant force 

To distinguish the dependences of the path probability on Lagrangian and 
Hamiltonian actions, it is necessary to random motion under conservative 
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Figure 1: The result of numerical simulation of stochastic motion with 10^ free 
particles. The left panel shows the sampled paths between the given 
points a and b. The right panel shows the path probability distribu- 
tion against the actions. The Lagrangian and Hamiltonian actions are 
equivalent as V{x) = for free particles. The straight line is a best fit 
of the points with a slope of - 6.7 x 10^'' J~^s~^ 

forces. The first force we studied is the constant force obtained from the 
potential V{x) = mgx. The regular motion is described by f{t) = —^gt"^, 
where the parameter g = lOm/s^. The results are shown in the right panel 
of Figure 2. Eq. (4) still holds with 7 f« 6.4 x 10^^ J~^s"^. There is no 
correlation between path probability and Hamiltonian action. 

3.3 Particles under harmonic force 

The potential of the harmonic force is V{x) = ^kx'^ giving a regular motion 
f{t) = Asin{ut), where A = Ix 10"^ m and u = {k/niY/'^ = 4.7 x 10^ s"^ 
The right panel of Figure 3 shows the path probability distribution against 
actions. As for constant force, the path probability distribution decreases 
exponentially with increasing Lagrangian action with a slope of the straight 
line — 6.6 x 10^^ J~^s~^. No correlation with the Hamiltonian action is 
found. 

3.4 Particles in cubic potential 

To our opinion, the above results with 3 potentials are sufficiently convincing 
for the claim that the path probability decreases exponentially with increas- 
ing the Lagrangian action instead of the Hamiltonian one. But by curiosity. 
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Figure 2: The result of numerical simulation of stochastic motion with 10 par- 
ticles subject to a constant force with potential V{x) = mgx. The left 
panel shows the different sampled paths between the given points a 
and b. The right panel shows the path probability distribution against 
the Lagrangian (circles) and Hamiltonian (stars) actions. The straight 
line is a best fit of the points. It implies an exponential dependence 
on the Lagrangian action with 7 6.4 x 10^^ J~^s~^ in Eq. (4). No 
correlation between the path probability and the Hamiltonian action is 
found. 
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Figure 3: The result of numerical simulation of stochastic motion with 10^ par- 
ticles subject to a harmonic force with V{x) = ^kx'^. The left panel 
shows the different sampled paths between the given points a and b. 
The left panel shows the different sampled paths between the given 
points a and b. The right panel shows the path probability distribu- 
tion against Lagrangian (circles) and Hamiltonian (stars) actions. The 
straight line is a best fit of the points whose slope gives 7 6.6 x 10^^ 
J-^s^i in Eq. (4). 
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Figure 4: The result of numerical simulation of stochastic motion with 10^ par- 
ticles subject to a cubic potential V{x) = ^Cx^. The left panel shows 
that the different sampled paths between the given points a and b. The 
right panel shows the path probability distribution against Lagrangian 
(circles) and Hamiltonian (stars) actions. The straight line is a best fit 
of the points whose slope gives 7 4.5 x 10^^ J~^s^^ for Eq. (4). 

we also tried two other higher order potentials. The first one is V{x) = ^Cx^ 
giving a regular motion f{t) = -^(§T^F (^0 = 3 x 10-^ C = 200). The 
path probability distributions against the two actions are shown in Figure 4. 
Eq. (4) holds for the Lagrangian action with the coefficient 7 ~ 4.5 x 10^^ 

3.5 Particles in quartic potential 

For the one-dimensional quatric oscillator [17, 18, 19], the potential has the 
form V{x) = \Cx'^, with a regular motion f{t) = AmSin{wt). Unlike the 
harmonic potential, the frequency w will depend on the amplitude which is 
= Ix IQ-^ m, giving u = ^ ^ (|g)^/M^ = 2 x 10^ s'^ [17], where 
T is the complete cycle period. The path probability distributions against 
the two actions are shown in Figure 5. The distribution Eq. (4) with the 
Lagrangian action is still confirmed with 7 ^ 7.8 x 10^^ J"^s~^. 

4 Discussion 

From the above figures, the exponential dependence of the path probability 
on the Lagrangian action is obvious. This correlation can be characterized 
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Figure 5: The result of numerical simulation of stochastic motion with 10^ par- 
ticles subject to a quartic potential. The left panel shows the differ- 
ent sampled paths between the given points a and b. The right panel 
shows the path probability distribution against Lagrangian (circles) and 
Hamiltonian (stars) actions. The straight line is a best fit of the points 
whose slope gives 7 7.8 x 10^'' J^^s"^ for Eq. (4). 



by a correlation function between A {A^ or Ah)) and — In P{A) given by 

^ EUiA- < Aj >)[-lnP{A,)+ < InPjA,) >] 

^J[EtliA- < A, >)^][Etii-lnP{A,)+ < lnP{A,) >)2] ' 

where < Ai > and < In P{Ai) > are the means of action A and — In P[A) 
respectively. c{A, — In P{A)) ^ 1 would indicate that A and —lnP{A) are 
linearly correlated. The results obtained from the numerical experiments are 
shown in Table 1. 

It can be concluded from the Table 1 that — In P{A) has a linear corre- 
lation with the Lagrangian action and no correlation with the Hamiltonian 
one. It should be indicated that all the numerical simulations have been done 
with different time scale of each step. No dependence of c{A) on time scale 
was observed. 

It can be expected that the result would be more precise with less un- 
certainty if we could use more particles and more time steps. More precise 
simulation will be possible with more powerful computer. 
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Table 1: Values of the correlation function c{A) between the logarithm of the 
path probability — InP(^) and the Lagrangian action Al in comparison 
with the Hamiltonian one Ah for the 5 potentials V{x) used in the 
numerical experiments. The linear correlation between — InP(^) and 
Al is obvious. 



V{x) c{Al) c{Ah) 

0.9865 0.9865 

mgx 0.9686 0.4206 

Ikx^ 0.9827 0.2975 

ICx^ 0.9162 0.2302 

\Cx^ 0.9397 0.5635 



5 Conclusion 

To summarize, we made numerical experiments of stochastic motion of non 
dissipative systems or weakly dissipative systems in order to calculate the 
path probability and to investigate its dependence the conventional mechan- 
ical quantities. The model of the simulation is small silica particles subject 
to conservative forces and Gaussian noises. The time scale of observation, or 
of each step of the simulation is larger than 10~^ s (simulations with smaller 
time scale of observation and non Gaussian noises will be described later). 
It is found that the path probability decreases exponentially with increasing 
action (Lagrangian one) of the paths, meaning that the most probable path 
is just the least action path of Hamiltonian/Lagrangian mechanics. This 
is a reasonable result of the model since, with diminishing noise, more and 
more paths shrink onto the bundle of least action paths. In the limiting case 
of vanishing noise, all paths will collapse on the least action path and the 
motion will recover the Hamiltonian/Lagrangian dynamics. 

From this results, it can be predicted that, for such kind of motions, the 
probability of occurrence of any path from a given point to any arbitrarily 
chosen point, within a given duration of motion, must decrease exponentially 
with increasing action. An interesting case to be studied is the probability 
of the paths from one point to another with fixed energy but arbitrary time 
duration such as in the very first formulation of least action principle given 
by Maupertuis[20]. 



10 



Finally we would like to mention that, this exponential path probability 
is one of the possible path probability distributions underlying a stochastic 
formalism of Hamiltonian mechanics. The fundamental points of this for- 
mulation have been laid down in [7, 8] with several theoretical consequences 
described in [21]. 
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